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We show in the Skyrme-Hartree-Fock approach that unambiguous correlations 
exist between observables of finite nuclei and nuclear matter properties. Using 
this correlation analysis to existing data on the neutron skin thickness of Sn 
isotopes, we find important constraints on the value -Esym(po) and density 
slope L of the nuclear symmetry energy at saturation density. Combining these 
constraints with those from recent analyses of isospin diffusion and double 
neutron/proton ratio in heavy ion collisions leads to a value of L = 58 ± 18 
MeV approximately independent of -Esym{po)- 

Keywords: Nuclear symmetry energy; Neutron skin; Skyrme Hartree-Fock. 



1. Introduction 

The nuclear symmetry energy i?sym(p) plays a crucial role in both nuclear 
physics and astrophysics.^'^ Although significant progress has been made 
in recent years in determining the density dependence of i?sym(/o),^ large 
uncertainties still exist even around the normal density po,^ and this has 
hindered us from understanding more precisely many important properties 
of neutron stars. ^ To constrain the symmetry energy with higher accuracy 
is thus of crucial importance. 
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Theoretically, it has been established^"^'^ that the neutron skin thick- 
ness Ar„p = {r^Y^"^ — (j'p)"'^^^ of heavy nuclei, given by the difference of 
their neutron and proton root-mean-squared radii, provides a good probe of 
Egyn-iip)- In particular, Ar„p has been found to correlate strongly with both 
-E'sym(po) and L in microscopic mean-field calculations. ^^^^ It is, however, 
difficult to extract an accurate value for L from comparing the calculated 
Ar„p of heavy nuclei with experimental data as it depends on several nu- 
clear interaction parameters in a highly correlated manncr^'^ and the calcu- 
lations have been usually carried out by varying the interaction parameters. 
A well-known example is the Skyrme-Hartree-Fock (SHF) approach using 
normally 9 interaction parameters and there are more than 120 sets of 
Skyrme interaction parameters in the literature. 

In the present talk, we report our recent work^^ on a new method to 
analyze the correlation between observables of finite nuclei and some macro- 
scopic properties of asymmetric nuclear matter. Instead of varying directly 
the 9 interaction parameters within the SHF, we express them explicitly 
in terms of 9 macroscopic observables that are either experimentally well 
constrained or empirically well known. Then, by varying individually these 
macroscopic observables within their known ranges, we can examine more 
transparently the correlation of Ar„p with each individual observable. In 
particular, we have demonstrated that important constraints on Esymipo) 
and L can be obtained with the application of this correlation analysis to 
existing data on the neutron skin thickness of Sn isotopes. 



2. The theoretical method 

In the standard SHF model, the 9 Skyrme interaction parameters, i.e., 
(7, to ~ ^3- xq ~ X3 can be expressed analytically in terms of 9 macroscopic 
quantities po, Eq{pq), the incompressibility Kg, the isoscalar effective mass 
, the isovector effective mass m* q, £^sym(/3o), L, the gradient coefficient 



Gs, and the symmetry-gradient coefficient Gy,^^'^^ i.e., 

(7 = 7-1, io = 4a/(3po), xo = 3(y-l)i?^°J„(po)/a-l/2, 
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where fcO = {Lbn'^p^f'^ E'°Upo) - ^^symCpo) ^ 
parameters C, I?, a, /3, 7, and ?/ are defined as^^ 
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As a reference for the correlation analyses below, we use the MSLO 
parameter set,^^ which is obtained by using the following empirical values 
for the macroscopic quantities: po — 0.16 fm"-^, Eo{po) = —16 MeV, Kq = 
230 McV, m* o = 0.8to, „ = 0.7m, E^yn^ipo) ^ 30 McV, and L 60 
MeV, Gy 5 MeV•fm^ and Gs = 132 MeV-fm^. And the spin-orbital 
coupling constant Wq = 133.3 MeV -fm^ is used to fit the neutron pi/2 ~ 
P3/2 splitting in ^^O. Using other Skyrme interactions obtained from fitting 
measured binding energies and charge rms radii of finite nuclei does not 
change our conclusion. 



3(7 - l)i?'.%(Po) 
2/3 



-,(2) 



3. Results 

To reveal clearly the dependence of Ar„p on each macroscopic quantity, we 
vary one quantity at a time while keeping all others at their default values 
in MSLO. Shown in Fig. 1 are the values of Ar„p for ^ospb, ^^°Sn and ^^Ca. 
Within the uncertain ranges for the macroscopic quantities considered here, 
the Ar„p of ^°^Pb and ^^°Sn exhibits a very strong correlation with L. 
However, it depends only moderately On 
the other hand, the Ar„j, of ^*Ca displays a much weaker dependence on 
both L and Egymipo)- Instead, it depends moderately on Gy and Wq. This 
explains the weaker Ar„p-iJsym(p) correlation observed for ^*Ca in previous 
SHF calculations using different interaction parameters. 

Experimentally, the Ar„p of heavy Sn isotopes has been systematically 
measured. ^^"^'^ As an illustration, we first show in the panel (a) of left 
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Fig. 1. (Color online) The neutron skin thickness Ar„p of Pb, ^^"Sn and ^*Ca from 
SHF with MSLO by varying individually L (a), Gy (b), Gs (c), Eo{po) (d), Saym(po) 
(e), Ko (f), m*_o (g), m; (h), po (i), and Wo (j). Taken from Ref." 



window in Fig. 2 the comparison of the available Sn Ar„p data with our 
calculated results using 20, 60 and 100 MeV, respectively, for the value of 
L and the default values for all other quantities in MSLO. It is seen that 
the value L = 60 MeV best describes the data. To be more precise, the 
evaluated from the difference between the theoretical and experimental 
Ar„p values is shown as a function of L in the panel (b) of left window in 
Fig. 2. The most reliable value of L is found to be L = 54 ± 13 MeV within 
a 2(T uncertainty. 

Since the value of Ar„p depends on both L and i?sym(po), we have 
carried out a two-dimensional analysis as shown by the grey band in 
the panel (c) of left window in Fig. 2. It is seen that increasing the value 
of Egymipo) systematically leads to smaller values of L. Furthermore, we 
have estimated the effects of nucleon effective mass by using m*_g = 0.7m 
and m* q = 0.6m as well as m* q = 0.9m and m* q = 0.8m, in accord with 
the empirical constraint m* q > m* qj^'^'^^ and the resulting constraints are 
shown by the dashed and dotted lines. As expected from the results shown 
in Fig. 1, effects of nucleon effective mass are small with the value of L 
shifting by only a few MeV for a given £^sym(po)- As one has also expected, 
effects of varying other macroscopic quantities are even smaller. 



January 16, 2013 0:58 



WSPC - Proceedings Trim Size: 9in x 6in ChenLW-NSI<in 



5 

A L (MeV) 



112 116 120 124 128 132 40 60 80 




26 28 30 32 34 



Fig. 2. (Color online) Left window: (a) The Ar„p data for Sn isotopes from different 
experimental methods and results from SHF calculation using MSLO with L = 20, 60 
and 100 MeV. (b) as a function of L. (c) Constraints on L and -Bsym(po) from the 
analysis of the Ar„p data on Sn isotopes (Grey band as well as dashed and dotted lines). 
Right window: Constraints on L and -E'sym(po) obtained in the present work (dashed 
lines) and that from Ref.^^ (solid lines). The shaded region represents their overlap. 
Taken from Ref.^" 

The above constraints on the L-Esym{po) correlation can be combined 
with those from recent analyses of isospin diffusion and double n/p ratio in 
heavy ion collisions at intermediate cnergics^^ to determine simultaneously 
the values of both L and i?sym(po)- Shown in the right window of Fig. 2 
are the two constraints in the Egy^^poj-L plane. Interestingly, these two 
constraints display opposite L-Egymipo) correlations. This allows us to ex- 
tract a value of L = 58 ± 18 McV approximately independent of the value 
of Egymipo)- This value of L is quite precise compared to existing estimates 
in the literature (See Ref.^ for a recent summary) although the constraint 
on Egymipo) is not improved. 

4. Summary 

We have proposed a new method to analyze the correlations between ob- 
servables of finite nuclei and some macroscopic properties of nuclear mat- 
ter, and demonstrated that the existing neutron skin data on Sn isotopes 
can give important constraints on the symmetry energy parameters L and 
Esymipo)- Combining these constraints with those from recent analyses of 
isospin diffusion and double n /p ratio in heavy ion collisions leads to a quite 
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accurate value of L = 58 ±18 MeV approximately independent of i?sym(po)- 
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